Draft version February 5, 2008 

Preprint typeset using I^T^X style cmulateapj v. 08/13/06 



THE SLOAN LENS ACS SURVEY. IV: THE MASS DENSITY PROFILE OF EARLY-TYPE GALAXIES OUT 

TO 100 EFFECTIVE RADII 

Raphael Gavazzi 1 , Tommaso Treu 1 , Jason D. Rhodes 2,3 , Leon V.E. Koopmans 4 , Adam S. Bolton"', Scott Burles 6 , 

Richard J. Massey 3 , and Leonidas A. Moustakas 2 

Draft version February 5, 2008 



o 
o 

(N 
u 

< 

OO 



(N 
> 

OO 

in 
o 
o 

Or 

6 



ABSTRACT 

We present a weak gravitational lensing analysis of 22 early-type (strong) lens galaxies, based on 
deep Hubble Space Telescope images obtained as part of the Sloan Lens ACS Survey. Using advanced 
techniques to control systematic uncertainties related to the variable point spread function and charge 
transfer efficiency of the Advanced Camera for Surveys (ACS), we show that weak lensing signal is 
detected out to the edge of the Wide Field Camera (< 300 h^ 1 kpc at the mean lens redshift z = 0.2). 
We analyze blank control fields from the COSMOS survey in the same manner, inferring that the 
residual systematic uncertainty in the tangential shear is less than 0.3%. A joint strong and weak 
lensing analysis shows that the average total mass density profile is consistent with isothermal (i.e. 
p oc r~ 2 ) over two decades in radius (3-300 ft^kpc, approximately 1-100 effective radii). This 
finding extends by over an order of magnitude in radius previous results, based on strong lensing 
and/or stellar dynamics, that luminous and dark component "conspire" to form an isothermal mass 
distribution. In order to disentangle the contributions of luminous and dark matter, we fit a two- 
component mass model (de Vaucouleurs + Navarro Frenk & White) to the weak and strong lensing 
constraints. It provides a good fit to the data with only two free parameters; i) the average stellar 
mass-to-light ratio M^jLy = 4.48±0.46 /iMq/L© (at z = 0.2), in agreement with that expected for an 
old stellar population; ii) the average virial mass-to-light ratio M v - 1T /Lv = 246+gy 1 /iM Q /L Q . Taking 
into account the scatter in the mass-luminosity relation, this latter result is in good agreement with 
semi-analytical models of massive galaxies formation. The dark matter fraction inside the sphere of 
radius the effective radius is found to be 27±4%. Our results are consistent with galaxy-galaxy lensing 
studies of early-type galaxies that are not strong lenses, in the region of overlap (30-300 /i _1 kpc). 
Thus, within the uncertainties, our results are representative of early-type galaxies in general. 
Subject headings: gravitational lensing - dark matter - galaxies : Ellipticals and lenticulars, cD - 
galaxies: structure 



1. INTRODUCTION 

It is now commonly accepted that cold non-baryonic 
dark matter dominates the dynamics of the Uni- 
verse. Whereas the so-called ACDM (cold dark mat- 
ter) paradigm has been remarkably successful at repro- 
ducing with high level of precision the properties of the 
universe on scales larger than Mpc (e.g. Spergel et al. 
2006; Tegmark et al. 2004; Seljak et al. 2005), the sit- 
uation at galactic and subgalactic scales is more uncer- 
tain. Dark-matter-only numerical simulations make very 
clear predictions. Dark matter halos have a characteris- 
tically "cuspy" radial profile (e.g. Navarro et al. 1997; 
Moore et al. 1998; Ghigna et al. 1998; Jing 2000; Stoehr 
et al. 2002; Navarro et al. 2004), are triaxial (Jing 2002; 
Kazantzidis et al. 2004; Hayashi et al. 2006), and have 
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abundant substructure (e.g. Moore et al. 1999; De Lucia 
et al. 2004; Gao et al. 2004; Taylor & Babul 2004). From 
an observational point of view, substantial effort has been 
devoted to comparing those predictions to observations 
with debated results in the case e.g. of low surface bright- 
ness galaxies (Salucci 2001; de Blok et al. 2003; Swaters 
et al. 2003; Gentile et al. 2004; Simon et al. 2005) or at 
galaxy cluster scales (e.g. Sand et al. 2004; Gavazzi 2005; 
Comerford et al. 2006). The main source of ambiguity 
in such comparisons is due to the effects of baryons. Al- 
though a minority in terms of total mass, baryons are dis- 
sipative and spatially more concentrated than the dark 
matter, playing a critical role at scales below tens of kilo- 
parsecs. Understanding baryonic physics and the inter- 
play between dark and luminous matter is necessary to 
understand how galaxies form and, ultimately, to test 
the cosmological model. From an observational point of 
view, measuring the relative spatial distribution of stars, 
gas, and dark matter is essential to provide clues to help 
understand the physical processes and hard numbers to 
perform quantitative tests of models. 

The dark halos of early-type (i.e. elliptical and lentic- 
ular) galaxies have been particularly hard to detect and 
study, due to the general lack of kinematic tracers, such 
as HI, at large radii. Studies of local galaxies based on 
stellar kinematics (Bertin et al. 1994; Gerhard et al. 2001; 
Mamon & Lokas 2005a, b; Cappellari et al. 2006), kine- 
matics of planetary nebulae (Romanowsky et al. 2003; 
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Arnaboldi et al. 2004; Merrett et al. 2006) and temper- 
ature profile of X-ray emitting plasma (Humphrey et al. 
2006) indicate that at least for the most massive sys- 
tems dark matter halos are generally present. The total 
mass density profile is found to be close to isothermal 
(ptot ex r ~ 2 ) °n scales out to a few effective radii. 

In the distant universe an additional mass tracer is pro- 
vided by gravitational lcnsing. At scales comparable to 
the effective radius, strong gravitational lensing makes it 
possible to detect and study the mass profile and shape 
of individual halos (Kochanek 1995) or statistically of a 
population of halos (Rusin et al. 2003). The combina- 
tion of strong lensing with stellar kinematics (Treu & 
Koopmans 2002; Koopmans & Treu 2002, 2003; Treu & 
Koopmans 2004; Koopmans et al. 2006) is particularly 
effective, and allows one to decompose the total mass 
distribution into a luminous and dark matter with good 
precision, yielding information on the internal structure 
of early-type galaxies all the way out to the most distant 
lenses known at z <~ 1. At larger scales, the surface mass 
density is too low to produce multiple images. However, 
the weak lensing signal can be detected statistically by 
stacking multiple galaxies and measuring the distortion 
of the background galaxies. The statistical nature of this 
measurement imposes a certain degree of spatial smooth- 
ing or binning, which in turn limits the angular resolution 
of weak lensing studies. 

In this paper we exploit deep ACS images of 22 gravi- 
tational lenses from the Sloan Lens A CS Survey (Bolton 
et al. 2006; Treu et al. 2006a; Koopmans et al. 2006, 
hereafter papers I, II and III, respectively) to perform a 
joint weak and strong lensing analysis. This allows us 
to bridge the gap between the two regimes and study 
the mass density profile of early-type galaxies across the 
entire range ~1-100 effective radii, disentangling the lu- 
minous and dark components. 

From a technical point of view, joint weak and strong 
lensing analysis has already been applied in the past to 
clusters of galaxies and galaxies in clusters (Natarajan & 
Kneib 1997; Geiger & Schneider 1998; Natarajan et al. 
2002; Kneib et al. 2003; Gavazzi et al. 2003; Bradac et al. 
2005b, a). However, there are important differences with 
respect to galaxy scales. First and foremost, clusters pro- 
duce a much stronger weak lensing signal and therefore 
it can be detected and studied for individual systems. 
Second, since clusters are spatially more extended than 
galaxies, relatively large smoothing scales can be adopted 
to average the signal over background galaxies. In con- 
trast, the signal of individual galaxies is too "weak" to be 
detected, so that stacking a number of lens galaxies is re- 
quired to reach a sufficient density of background sources. 
For this purpose, previous studies have typically relied on 
very large sample of galaxies (Brainerd et al. 1996; Grif- 
fiths et al. 1996; Wilson et al. 2001; Guzik & Seljak 2002; 
Hoekstra et al. 2004, 2005; Klcinheinrich et al. 2006) or 
recently in the SDSS survey (Sheldon et al. 2004; Man- 
delbaum et al. 2006, hereafter S04 and M06 respectively). 
As we will show in the rest of the paper, the high den- 
sity of useful background galaxies afforded by deep ACS 
exposures (~72 per square arcmin) allows us to achieve 
a robust detection of the weak lensing signal with only 
22 galaxies, and to study the mass density profile with 
unprecedented radial resolution. 

The paper is organized as follows. After briefly summa- 



rizing the gravitational lensing formalism and notation in 
§ 2, we discuss the sample, data reduction and analysis, 
and the main observational properties of the lens galaxies 
in § 3. § 4 details the shear measurement, with a particu- 
lar emphasis on the precision correction for instrumental 
systematic effects and on tests of residual systematics by 
means of a parallel analysis of blank fields. This section 
also presents the mean radial shear profile around SLACS 
strong lenses and high resolution two-dimensional mass 
reconstructions. We combine strong and weak gravita- 
tional lensing constraints in § 5 to model the radial pro- 
files lenses and disentangle the stellar and dark Matter 
components. We discuss our results in § 6 and give a 
brief summary in § 7. 

Throughout this paper we assume the concordance cos- 
mological background with Hq = 100 h kms" 1 Mpc -1 , 
Sl m = 0.3 and f^A = 0.7. All magnitudes are expressed 
in the AB system. 

2. BASIC LENSING EQUATIONS 

In this section we briefly summarize the necessary 
background of gravitational lensing and especially the 
weak lensing regime which concerns the present analy- 
sis. The main purpose of this section is to define no- 
tations. We refer the reader to the reviews of Mcllicr 
(1999), Bartelmann & Schneider (2001) and Schneider 
(2006) for more detailed accounts. 

The fundamental quantity for gravitational lensing is 
the lens potential ip(9) at angular position 9 which is 

related to the surface mass density £(#) projected onto 
the lens plane through 

m = f dV£(0')m|0-0'| , (1) 

where D Q \, D os and D\ s are angular distances to the lens, 
to the source and between the lens and the source respec- 
tively. The deflection angle a. = Vtp relates a point in 
the source plane (3 to its image(s) in the image plane 9 
through the lens equation (3 = 9—a{9). The local relation 
between (3 and 9 is the Jacobian matrix ay = d(3i/d9j 

('"^ i-7; 71 ) • ( 2 ) 

The convergence k{9) = E(#)/£ cr i t is directly related 
to the surface mass density via the critical density 



and satisfies the Poisson equation 

Atp = tp.ii + ^,22 = 2k . (4) 

The 2-component shear is 7 = 71 + 272 = 5(^,11 — 
V',22) + iip, 12 in complex notation. An elliptical object 
in the image plane is characterized by its complex cllip- 
ticity e. In the weak lensing regime, the source intrinsic 
cllipticity e s and e are simply related by e = e s + 7. 

It is convenient to express the shear into a tangential 
and a curl term 7 = jt+ilx such that j t = — H { 7 e~ 24V } 
and 7 X = — X {je~ 2lip } with ip the polar angle. For 
a circularly symmetric lens, 7 X vanishes whereas 74 at 
radius r can be written as the difference between the 
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mean convergence within that radius k(< r) and the local 
convergence at the same radius n(r) : 

7t = «(< r) - n(r) . (5) 

In equations 1 to 5 we can isolate a geometric term which 
linearly scales the lensing quantities k, ip, and 7 and 
only depends on the distance ratio D\ s /D os . We thus 
can write ip — w(z\, Zg^oo (and analogously for k and 
7) with w(zi, z s ) = D\ s / D os Q(z s — z\), where Q(x) is the 
Heaviside step function. If sources are not confined in a 
thin plane, we account for the distribution in rcdshift by 
defining an ensemble average distance factor w(zi) such 
that: 



w(zi) = (w(zi,z s )) Ze = 




/ dz s n(z s ) 
Jo 

3. THE DATA 
3.1. Lens sample 

In this paper we focus on a subsample of 22 lens early- 
type galaxies taken from the SLACS Survey (paper I). 
The subsample is defined as all the confirmed lenses 
for which deep, 1-orbit long, ACS image through filter 
F814W were available as of the cutoff date for this pa- 
per, 2006 October 15. The parent sample is spectroscop- 
ically identified from the SDSS database and confirmed 
by ACS imaging, as described in paper I (see also Bolton 
et al. 2004). Ten lenses are in common with the sam- 
ple previously analyzed in papers II and III, while the 
remaining 12 lenses were not analyzed in papers II and 
III. Images of the 12 new lenses are shown in Figure 1. 
Table 1 summarizes the most relevant properties of the 
22 lenses. More details on the new lenses and on the 
ongoing programs will be presented elsewhere. 

The SDSS aperture velocity dispersions arc in the 
range 196 < a v < 333 km s -1 and the mean square veloc- 
ity dispersion is (c^} 1//2 — 248 kms -1 . The lens galax- 
ies have a mean redshift (z\ = 0.22). Paper II showed 
that SLACS strong lenses fall on the same Fundamen- 
tal Plane of non-lens early-type galaxies (see also Bolton 
et al. 2007, in prep). This demonstrates that - within 
our measurement errors - lensing galaxies have normal 
internal dynamical properties at small scales. One of 
the goals of this paper is to combine strong with weak 
lensing to check whether the outer regions of the SLACS 
lenses behave in a peculiar way as compared to non-lens 
early- type galaxies. 

3.2. HST observations & data reduction 

The sample was observed with the Advanced Camera 
for Surveys on board HST between November 2005 and 
October 15 2006, as part of HST programs 10494 (PI: 
Koopmans) and 10886 (PI: Bolton). One-orbit exposures 
were obtained through filter F814W (hereafter /) with 
the Wide Field Camera centering the lens on the WFC1 
aperture, i.e. in the center of the second CCD. Hence the 
observations cover a region as far as 3' around the lenses. 
Four sub-exposures were obtained with a semi-integer 
pixel offset (acs-wf c-dither-box) to ensure proper cos- 
mic ray removal and sampling of the point spread func- 
tion. For the lenses in program 10494 additional one or- 
bit exposures with ACS through filter F555W and with 



the NICMOS NIC2 camera through filter F160W are also 
available. In this paper, the additional F555W exposure 
is used to check satellite/foreground contamination to 
the weak-lensing catalog. A full detailed analysis of the 
multi-color images will be presented elsewhere. 

Since the goal of this paper is detecting the weak lens- 
ing signal produced by the SLACS strong lenses, we op- 
timize our reduction according to the prescriptions of 
Rhodes et al. (2007, hereafter R07). For each target, we 
used multidrizzle (Koekemoer et al. 2002) to combine 
the four subexposures, using a final pixel size of 0'.'03 and 
a Gaussian interpolation kernel. 

One important difference between this study and that 
of R07, is that ours are pointed observations. Thus in- 
strumental effects could play a different role than for 
weak-lensing studies of objects at random positions on 
the detector, possibly introducing systematic errors. In 
order to determine the amplitude of systematic errors in 
our weak lensing analysis, we carried out a perfectly anal- 
ogous analysis of 100 / images of the COSMOS survey 7 
with identical depth. As detailed below, this allows us to 
infer the mean and field-to-field variance of instrumental 
biases, showing that they are negligible for our purposes. 

3.3. Surface photometry and lens models 

Surface photometry of the lens galaxies was obtained 
by fitting de Vaucouleurs profiles after carefully masking 
the lensed structures (rings) and any neighboring bright 
satellites. The two-dimensional parametric fit was car- 
ried out using galfit (Peng et al. 2002). We checked 
that our results are consistent with those of paper II for 
the 10 lenses previously observed with shallower HST 
snapshot imaging (see corrected table 2 of paper II in 
Treu et al. 2006b). 

We determined absolute V band magnitudes of the 
lenses taking into account filter transformations and 
galactic extinction according to the Schlcgel et al. (1998) 
dust maps. Furthermore, in order to homogenize the 
sample, we passively evolved all V band magnitudes to 
a fiducial rcdshift z = 0.2, using the relation (Treu et al. 
(2001) and paper II): 

dlog^ 

— ^- ~ -0.40 ± 0.05 , (7) 

dz 

which is well suited for the massive early-type galaxies 
we are considering here. We note that the correction is 
of order a few hundredths dex, and adopting a differ- 
ent passive evolution would not alter our results in any 
significant way. Thus the V-band luminosities listed in 
Table 1 are z = 0.2 V band luminosities and can be con- 
sidered as fair proxies for the lens stellar mass up to an 
average stellar mass-to-light ratio Ty = M*/Ly. 

We measured Einstein radii in full analogy to paper 
III, that is, we parameterized the lens potential with a 
Singular Isothermal Ellipse profile and reconstructed the 
unlcnsed source surface brightness non-paramctrically to 
match the observed Einstein ring features. Typical un- 
certainties on the recovered values of i?Ein are <~ 0.05 
with small variation between lenses. Again, we checked 
that the present modeling provides consistent results 
with respect to those in paper III. A more detailed de- 

7 http : //www . astro . caltech . edu/~ cosmos/ 
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Fig. 1. — F814W gallery of the 12 new SLACS lenses. The postage stamps are 7.2 arcseconds on a side. From top 
to bottom: SDSSJ0029-0055, SDSSJ0157-0056, SDSSJ0252+0039, SDSSJ0330-0020, SDSSJ0728+3835, SDSSJ0808+4706, 
SDSSJ0903+4116, SDSSJ1023+4230, SDSSJ1103+5322, SDSSJ1205+4910, SDSSJ2238-0754, and SDSSJ2341+0000. For each 
lens we show the original image (left) and after subtraction of a model for the lens surface brightness (right). 



scription of strong lcnsing modeling of the 12 new lenses 
will be given in forthcoming papers. 

4. SHEAR ANALYSIS 

4.1. Background sources selection 

The detection of background sources was done 
with imcat 8 and cross-correlated with the SExtractor 
(Bertin & Arnouts 1996) source catalog to remove spu- 
rious detections. To limit screening by the foreground 
main lens we subtracted its surface brightness profile be- 
fore source detection with SExtractor and imcat. After 
identifying stars in the magnitude size diagram in a stan- 
dard manner, we applied the following cuts to select ob- 
jects for which shapes could be reliably measured. First, 
we restricted the analysis to galaxies brighter than I < 26 
although the galaxy sample is complete down to / ~ 27.5, 
based on the number counts. This removes faint and 
small objects with poorly known redshift distribution. 
Second, we applied a bright / > 20 cut to the source 
sample, to minimize foreground contamination. Third, 
we discarded objects with an half-light radius r^ < 0.09" 
(for comparison the / PSF has r/j ~ 0.06"). Fourth, pairs 
of galaxies with small angular separation (< 0.5") were 
discarded since their shape cannot be reliably measured. 
After these cuts, we achieve a number density of useful 
background sources n^g — 72 arcmin~ 2 . 

8 http://www. if a.hawaii . edu/~kaiser/imcat/ 



The redshift distribution of sources is taken from the 
recently published COSMOS sample of faint galaxies de- 
tected in the ACS/F814W band (Leauthaud et al. 2007). 
Their analysis exploits photometric redshifts, to rerive 
redshift sources distribution down to magnitude / ~ 26 
for a sizeable sample selected at HST resolution. The 
redshift distribution of sources having I< 26 is well rep- 
resented by the following expression 



dn(z s ) 
dz s 



1 



z T(a) 



e-^^izs/zoT- 1 



(8) 



with z = 0.345 and a = 3.89. For this particular red- 
shift distribution, values of %l = {D\ s / D os ) Zs are re- 
ported in Table 1. This redshift distribution represent 
a clear improvement of our analysis over previous esti- 
mates based on ground based surveys, as the redshift 
distribution of faint sources depends not only on magni- 
tude but also on object size. The relatively low redshifts 
of SLACS lenses and the rapid saturation of w(z\,z s ) 
with increasing source redshift helps reduce the sensitiv- 
ity of our results to residuals errors on photometric red- 
shifts. Taking into account current errors on dn(z s )/dz s 
reported by Leauthaud et al. (2007), the overall calibra- 
tion for our sample is accurate to ~ 7%. In the rest of 
the paper we will show that this uncertainty is negligible 
for our purposes. 

A potential additional concern is residual contamina- 
tion by satellites galaxies that are spatially correlated 
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TABLE 1 

Lens sample observed with deep ACS F814W imaging. 



Name 




Prog. ID 


Exp. time 


21 


0"ap 




Refi 


raj 


M v 


z B 


WSL 










(sec) 




(kms -1 ) 


(arcsec) 


(arcsec) 




+5 log h 








SDSSJ002907.8- 


005550 


10886 


2088 


0.227 


228 ± 18 


0.82 


1.48 


17.36 


-21.53 


0.931 


0.706 


0.721 


SDSSJ015758.9-005626 


10886 


2088 


0.513 


295 ± 47 


0.72 


0.93 


18.76 


-22.16 


0.924 


0.380 


0.441 


SDSSJ021652.5-081345 


10494 


2232 


0.332 


333 ± 23 


1.15 


2.79 


16.88 


-22.95 


0.523 


0.333 


0.608 


SDSSJ025245.2+003958 


10886 


2088 


0.280 


164 ± 12 


0.98 


1.69 


17.84 


-21.67 


0.982 


0.656 


0.662 


SDSSJ033012.1-002052 


10886 


2088 


0.351 


212 ± 21 


1.06 


1.17 


18.20 


-21.86 


1.107 


0.613 


0.589 


SDSSJ072805.0- 


f 383526 


10886 


2116 


0.206 


214 ± 11 


1.25 


1.33 


16.95 


-21.80 


0.688 


0.660 


0.745 


SDSSJ080858.8- 


f470639 


10886 


2140 


0.219 


236 ± 11 


1.23 


1.65 


17.10 


-21.77 


1.025 


0.735 


0.730 


SDSSJ090315.2- 


f411609 


10886 


2128 


0.430 


223 ± 27 


1.13 


1.28 


18.19 


-22.25 


1.065 


0.521 


0.512 


SDSSJ091205.3- 


f 002901 


10494 


1668 


0.164 


326 ± 12 


1.61 


5.50 


15.20 


-22.95 


0.324 


0.472 


0.794 


SDSSJ095944.1- 


f 041017 


10494 


2224 


0.126 


197 ± 13 


1.00 


1.99 


16.61 


-20.94 


0.535 


0.738 


0.840 


SDSSJ102332.3- 


f423002 


10886 


2128 


0.191 


242 ± 15 


1.30 


1.40 


19.93 


-21.56 


0.696 


0.686 


0.762 


SDSSJ110308.2- 


f532228 


10886 


2156 


0.158 


196 ± 12 


0.84 


3.22 


16.02 


-22.02 


0.735 


0.749 


0.801 


SDSSJ120540.4+491029 


10494 


2388 


0.215 


280 ± 13 


1.04 


1.92 


16.76 


-22.00 


0.481 


0.521 


0.735 


SDSSJ125028.3+052349 


10494 


2232 


0.232 


252 ± 14 


1.15 


1.64 


16.78 


-22.17 


0.795 


0.662 


0.716 


SDSSJ140228.1+632133 


10494 


2520 


0.205 


267 ± 17 


1.39 


2.29 


16.44 


-22.20 


0.481 


0.543 


0.747 


SDSSJ142015.9+601915 


10494 


2520 


0.063 


205 ± 04 


1.04 


2.49 


14.93 


-21.04 


0.535 


0.867 


0.919 


SDSSJ162746.5-005358 


10494 


2224 


0.208 


290 ± 14 


1.21 


2.47 


16.79 


-22.06 


0.524 


0.570 


0.743 


SDSSJ163028.2+452036 


10494 


2388 


0.248 


276 ± 16 


1.81 


2.01 


16.76 


-22.31 


0.793 


0.639 


0.698 


SDSSJ223840.2-075456 


10494 


2232 


0.137 


198 ± 11 


1.20 


2.33 


16.20 


-21.58 


0.713 


0.776 


0.827 


SDSSJ230053.2- 


f 002238 


10494 


2224 


0.228 


279 ± 17 


1.25 


2.22 


16.91 


-22.06 


0.463 


0.476 


0.719 


SDSSJ230321.7- 


H42218 


10494 


2240 


0.155 


255 ± 16 


1.64 


3.73 


15.97 


-22.40 


0.517 


0.670 


0.805 


SDSSJ234111.6- 


f 000019 


10886 


2088 


0.186 


207 ± 13 


1.28 


3.20 


16.30 


-22.14 


0.807 


0.729 


0.768 



Apparent I band magnitudes arc not corrected for Galactic extinction. Absolute magnitudes are K-corrected, extinction corrected, and 
corrected to the sample mean redshift z = 0.2 for luminosity evolution using logLy,0.2 = logLy,z + 0.40 * (z — 0.2). Combining 
measurement errors and uncertainties in various photometric corrections yields a typical error in apparent (resp. absolute) magnitudes 
±0.02 (resp. ±0.04) mag. Relative uncertainties in i?Ein are about 5% and ~ 10% for R c s- Since systems are elliptical, both i?Ein ami 
i? e ff are expressed relative to the geometric mean intermediate radius. i«sl is the lensing distance ratio for the strong lensing event source 
redshift whereas «J\vl is the same distance ratio averaged over the redshift distribution of background sources used for weak lensing. 



with the main lens galaxy and thus could dilute the weak 
lensing signal. Furthermore, we expect the number of 
satellites to depend on the distance from the lens cen- 
ter, and this could potentially affect our inferred shear 
profile. 

As a first check, we applied a color cut to the back- 
ground catalog of the 10 SLACS fields for which F555W 
filter imaging is available. We measured the shear sig- 
nal for galaxies redder than the lenses (i.e. — 
F814W > 1.5), expected to be at higher redshift (see 
e.g. Broadhurst et al. 2005; Limousin et al. 2006, for 
similar color selections). The signal-to-noise on the re- 
covered shear profile for this tiny subsample of sources 
turned out to be too small for this test to be conclusive. 
This test will be more powerful when the full multi-color 
dataset will be available at the end of the survey. 

Therefore we turned to comparisons with the weak 
lensing SDSS analysis of S04 who found that about 10% 
of r < 22 sources are correlated to the lens at scales 
~ 30kpc. Our ACS catalogs are 4 magnitudes deeper 
than SDSS catalogs. The number density of background 
sources is much higher N hg (I < 26)/N hg (I < 22) ~ 46 
but the number of satellites should also increase. As- 
suming a typical luminosity function with slope a — 
— 1, we can extrapolate our counts and predict that 
N S3jt (I < 26)/iV sat (7 < 22) < 2. Therefore, at ~ 30kpc 
from the lens center, the contamination must be at most 
Asat/Nbg — 10 x 2/46 ~ 0.5%. Similarly, at smaller scale 
r <~ 3 kpc, we can extrapolate S04 results to predict that 
the contamination ratio increases by at most a factor 10, 
yielding N sat /N^ g ~ 5%. This ratio, as we will see below, 
is much smaller than present error bars (> 30% per bin) 
so we conclude that contamination by satellites cannot be 



a significant source of error. This finding is supported by 
the excellent agreement between strong and weak lensing 
measurements at small scales (see below). 

4.2. Instrumental distortions 

Before using the shape of background galaxies as a 
tracer of the shear field, we need to correct several in- 
strumental effects. Because every lens galaxy is approx- 
imately 9 at the same location in the detector frame (in 
the middle of CCD2) we need to carefully assess and cor- 
rect any instrumental source of systematic polarization of 
galaxies which may bias the measured shear profile. To 
correct for the smearing of galaxy shape by the Point 
Spread Function (PSF), we use the well known KSB 
method (Kaiser et al. 1995) which has proved being a 
reliable method down to cosmic shear requirements (Hcy- 
mans et al. 2006b). The implementation we are using is 
similar to that of Gavazzi & Soucail (2006) but is tuned 
for the specific space-based conditions (see e.g. Hoek- 
stra et al. 1998) by adaptively matching the radial size 
of the weight function applied to stars to that of the 
galaxies that are being PSF-corrected. Some important 
changes inspired by R07 are detailed in the following (see 
also Schrabback et al. 2006, for further discussion of the 
techniques required to extract weak lensing signal from 
ACS images). 

4.2.1. Focus & Point Spread Function smearing 

The shape of galaxies must be corrected for the smear- 
ing by the Point Spread Function of the ACS camera 

9 Typically within a few pixels due to absolute pointing uncer- 
tainties. The stack is of course aligned on the measured center of 
each galaxy. 
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which circularizes objects and/or imprints systematic 
distortion patterns. The PSF from space-based images 
is expected to be more stable as compared to ground- 
based data which suffer from time varying atmospheric 
seeing conditions. However, R07 showed that the ACS 
PSF varies dramatically as a function of time, essentially 
because of focus oscillations due to thermal breathing. 
The peculiar off-axis position of ACS enhances any fo- 
cus variability and PSF anisotropy is difficult to control. 
Unfortunately, we cannot map PSF variations across the 
field from the data themselves, since not enough stars 
are observed in each exposure. Therefore, following R07, 
we compare the few available stars to mock PSFs built 
with TinyTim (Krist & Hook 2004) and modified as de- 
scribed in R07 as a function of focus and detetmine the 
best fitting focus. The distribution of offsets ^f OC us with 
respect to the nominal ACS focal plane is well consistent 
with R07 i.e. <5f ocus ~ —3.4 ± 0.8 fim. As a consistency 
check we apply the same procedure to the blank fields 
from COSMOS. 
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Fig. 2. — Ellipticity patterns measured onto stars in all the 
COSMOS (left) and SLACS (right) fields. Upper (black) pan- 
els show stellar ellipticities before applying PSF correction 
scheme whereas the lower (red) panels show that the main 
distortion patterns are significantly suppressed after correc- 
tion. Each panel corresponds to the full ACS field of view 
(i.e. 200" on a side) 

Fig. 2 shows the ellipticity of stars before and after 
our PSF correction scheme for the 100 COSMOS fields 
and our 22 SLACS fields. Averaged over the COSMOS 
fields we measured a mean complex ellipticity (e* )Uncor ) = 
(-0.0037 + 0.0072i) ± 0.0004 before PSF anisotropy cor- 
rection and (e*, cor ) = (0.0000 - 0.0043i) ± 0.0003 af- 
ter correction. Similarly, in SLACS fields we obtained 
(e*,uncor) = (-0.0025 + 0.0096i) ± 0.0006 and (e*, cor ) = 
(0.0011-0.0023i)±0.0005. In both datasets the scatter of 
corrected stellar ellipticities about this mean is isotropic 
and a e „ cor ~ 0.0087. We conclude that the correction re- 
duces the mean PSF anisotropy by a factor 2 — 4 although 
some residuals are still present at the ~ 0.003 level. We 



will show in § 4.2.3 that the residual uncertainties are 
negligible for our purpose. 

In addition to anisotropic distortions, the convolution 
with the PSF also produces isotropic smearing making 
objects appear rounder. This effect is much smaller on 
HST images than from the ground but it must be taken 
into account for small objects with a typical size compa- 
rable to that of the PSF. Our initial size cut 77, > 0.09" 
guarantees that such isotropic smearing is kept at a low 
level and the KSB method can perform an accurate cor- 
rection (Kaiser et al. 1995; Hoekstra et al. 1998). Our 
implementation of KSB builds on the proposed improve- 
ments suggested by the STEP1 and STEP2 results (Hey- 
mans et al. 2006b; Massey et al. 2007). These papers in- 
dicate that, in general, KSB methods can achieve ~10% 
relative shear calibration biases or smaller. Since this is 
smaller than our statistical errors, it is sufficient to adopt 
for the present paper a conservative 10% systematic un- 
certainty in our shear calibration (m STEP parameter). 
In a future paper we plan to take advantage of the future 
spaceSTEP simulation set to get a more accurate esti- 
mate of the uncertainty on the shear calibration. This 
will be necessary given the gain in sensitivity expected 
when the deep SLACS follow-up will be complete. In 
addition, we demonstrate in the next section that no 
significant residual additive term (c STEP parameter) 
is observed in either the SLACS data or the COSMOS 
control fields. 

4.2.2. Charge Transfer Efficiency 

Another source of systematic distortion is the degrada- 
tion of Charge Transfer Efficiency (CTE) on ACS CCDs. 
Charges are delayed by defects in the readout direction 
(i.e. the y axis, charges going from the gap between 
the CCD chips toward the field boundaries), imprinting 
a tail of electrons behind each object which modifies its 
shape, thus producing a spurious negative e\ component. 
Since the strength of CTE-induced distortions depends 
on the signal-to-noise (snr) of the objects (the fainter the 
source the stronger the effect), we cannot measure dis- 
tortions from stars and correct faint galaxies accordingly. 
This effect must be quantified and corrected with galaxies 
themselves but one must be able to distinguish between 
the physical signal and the CTE distortions. To this aim 
we use the blank COSMOS fields to make sure that our 
CTE correction scheme will efficiently remove CTE dis- 
tortions while leaving the actual shear signal unchanged. 
In practice, we use the empirical recipe proposed by R07 
in which an e^cTE component, function of y pixel co- 
ordinate, snr, observation time (since CCD degradation 
increases with time) is subtracted for each object. Here 
we adapt the expression from Eq. (10) of R07 



61 cte = -3.6 x 10" 



(MJD - 52333) ( 



snr 



0.18' 



-0.9 



-0.1 



(9) 



to the imcat definition of snr whereas R07 use 
SExtractor. Note that the small dependence of e^cTE 
on size is somewhat degenerate with the way snr is de- 
fined and may not be always necessary (like in R07). 
Note also that snr is calculated by imcat without taking 
into account noise correlation caused by multidrizzling 
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on oversampled pixels (0.03" pixel size instead of the na- 
tive 0.05" value). Therefore this expression may not be 
directly applicable with different multidrizzle settings. 
MJD is the Modified Julian Date of observation, rh is the 
half-light radius and y' is the normalized y pixel coordi- 
nate 10 . 

Fig. 3 illustrates the effect of CTE on e\ elliptic- 
ity components for COSMOS and SLACS fields (respec- 
tively left and right panels). We show the e\ component 
of galaxy ellipticity as a function of the y' frame coordi- 
nate for sources brighter than 1= 27 (i.e. well beyond our 
magnitude cut for selecting suitable sources). The upper 
panels show the mean e\ before and after CTE(+PSF) 
correction. For COSMOS and SLACS we see a similar 
tendency for vertical stretching of galaxies in the mid- 
dle of the frame. The empirical CTE distortion fitting 
formula (9) provides a good correction. Although sta- 
tistical errors are larger in the SLACS images (5 times 
smaller sample) we see a modulation of the corrected 
ei component as a function of y' that is not observed 
in the corrected COSMOS images. As we shall see be- 
low this is the signature of the signal we are interested 
in and it should not be erased by the CTE correction 
scheme (in order to compare this residual with the ex- 
pected shear, we overlay the [e\(y')) shear signal from 
an isothermal sphere with Einstein Radius i?Ein = 1.2"). 
The lower panels (left and right) show how CTE distor- 
tions depend on signal to noise ratio. We have split the 
galaxy sample into three magnitude quantiles, the bright 
objects being less distorted. This is well accounted for 
by the snr (and size) dependence in Eq. (9). The over- 
all amplitude of CTE distortion is approximately double 
for SLACS images because of increasing CCD degrada- 
tion with time. This is also well captured by Eq. (9). 
The median observation date of the 100 COSMOS expo- 
sures we are considering is MJD cosmos — 53141 and for 
SLACS it is WTDsUcs = 53972. 

4.2.3. Checks on residuals 

To test the quality of the instrumental systematics cor- 
rection, we plot in Fig. 4 the radial profile of both the 
tangential e t and curl e x components of the complex el- 
lipticity in COSMOS and SLACS fields (left and right 
respectively). The center is set on the lens for SLACS 
images and at the same location in the detector frame in 
COSMOS. If we first focus on these latter images, we see 
no statistically significant residual et nor e x component, 
thus showing that we are free from PSF (seen in stars 
and galaxies) or CTE (seen in galaxies only) systemat- 
ics. As well around SLACS lenses, stars do not carry 
any significant residual e t or e x signal. Therefore we can 
safely assume that our systematics correction scheme is 
accurate enough for the present analysis 11 . Galaxies in 
SLACS fields clearly carry a strong e t signal (note the 
first data point well outside the plotting window) whereas 
no statistically significant e x component is observed as 
expected for a gravitational lensing origin of this shear 
signal. 

10 y' = (rcsp. 1) at the bottom (resp. top) edge of the ACS 
field of view. 

11 At the end of the SLACS deep survey, we expect ~80-100 
lenses. Since systematics are already below statistical errors in 
COSMOS (100 fields), our treatment is satisfactory for the final 
sample. 



4.3. Other sources of error 

A final additional potential source of systematic uncer- 
tainy is the effect of the lens galaxy surface brightness on 
the ellipticity of background sources at small projected 
radii. To estimate this effect, we detected and measured 
object shapes before and after subtraction of the lens 
surface brightness profile using galf it and b-spline tech- 
niques developed for the strong lensing analysis (paper 
I). These two method yield similar results for the purpose 
of weak lensing. It turns out that the lens-subtraction 
process changes measured ellipticities by at most 5% in 
an incoherent way. Therefore we do not consider further 
the effect of lens surface brightness as a relevant potential 
source of systematic. 

4.4. Two-dimensional mass reconstruction and shear 

profile 

In the upper panel of Fig. 5 a mass reconstruction 
(convergence k map) around the stacked lenses using the 
Kaiser & Squires (1993) method is shown. There is only 
one significant convergence peak at the position of the 
main lens. Note that the Gaussian smoothing scale of the 
convergence maps is 8 arcsec. This extraordinary high 
spatial resolution is made possible by the high density of 
background sources. The lower panel shows the imagi- 
nary part of the reconstructed mass map (obtained after 
rotating background galaxies by 45°). This is consistent 
with a pure noise map and illustrates the amplitude of 
the noise. 

We now analyze the radial shear profile achieved by 
stacking the lens galaxies. Since the images are taken 
at a random orientation with respect to the lens major 
axis we can safely assume circular symmetry in the anal- 
ysis. We convert the shear 7 into the physical quantity 
AY,(R) = S crit 7(i?) = S(< R) - For a given 

lens redshift z\ we also define the average critical density 
s crit = -£g dJe^y An estimator for AS at a given 
radius is 

j = l ^critj L^i=\ ^ej , s 

critj L^di—\ u e,i 

where N s j is the number of sources in the radial bin 
around lens j and o~ e ^ is the uncertainty assigned to the 
tangential ellipticity estimate e*^ (see Gavazzi & Soucail 
2006, for details on this weighting scheme). With this 
definition, AS is directly comparable to other SDSS weak 
galaxy-galaxy lensing analyses (e.g. S04, M06). 

Measured AS values around SLACS lenses arc re- 
ported in Tab. 2 and shown in Fig 6. The detection 
significance is derived as follows. The x 2 of the data 
with respect to the zero shear hypothesis is 47.8 for 9 de- 
grees of freedom. The probability of finding a higher % 2 
is 3 x 10~ 7 , thus the non-detection hypothesis is rejected 
at the 99.99997% level. For a Gaussian distribution this 
is equivalent to 5a. 

To compare with previous studies, we consider the 
measurement from S04 for their subsample of massive 
a v > 186 km s -1 . The mean square velocity of their sam- 
ple is (ct 2 ) 1 / 2 ~ 225 km s~ 1 . In order to compare with our 
points we need to correct for the different velocity disper- 
sion. Assuming an isothermal profile, the shear scales as 
the velocity dispersion squared, so that we need to scale 
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Fig. 3. — Vertical spurious stretching (negative ei ellipticity component) as a function of vertical pixel coordinate due to 
CTE degradation. Measurements for COSMOS and SLACS fields are shown in the left and right panels respectively. Upper 
panels: average stretching before (solid black) and after (dashed red) applying our empirical correction term. Distortions are 
maximum far from the readout (gap between CCDs at y' ~ 0.5). They are stronger in SLACS images because the number of 
defects in CCDs increases with time. The corrected residual signal should not perfectly "flat" in SLACS fields because of the 
presence of the true weak lensing signal we aim at detecting. The amplitude of the mean ei shear signal expected from an 
isothermal sphere with Einstein Radius R^m = 1-2" is shown for comparison (solid blue). The lower left and lower right 
panels show that distortions are stronger for fainter objects as illustrated by splitting galaxies into three magnitude bins (red, 
green, blue curves for bright, intermediate and faint objects). Comparing ei ellipticities before (left) and after (right) applying 
our correction scheme, we see that equation (9) corrects this spurious signal at all magnitudes. 



their points up by (248/225) 2 ~ 1.21 for a proper com- 
parison. After this correction, the agreement is excellent 
in the radial range 60 < R < 200/i _1 kpc of overlap as 
shown in Figure 6. We also check that our results are 
in agreement with the AE profile of the sm7 (early-type) 
stellar mass bin of M06. 

5. JOINT STRONG & WEAK LENSING MODELING 

In this section we take advantage of the availability 
of both strong and weak lensing constraints to investi- 
gate the mass profile of SLACS lenses from a fraction 
of the effective radius to 100 effective radii (3 < R < 
300/i- x kpc). 

5.1. SIS consistency check 

Before considering more sophisticated models for the 
density profile, we first check if the singular isothermal 
density profile favored in the inner parts of galaxies by 
strong lensing alone (e.g. Rusin et al. 2003), and by 
strong lensing and stellar dynamics (Treu & Koopmans 
2004, paper III), is consistent with our weak lensing mea- 
surements. 

5.1.1. Consistency with strong lensing 

For a singular isothermal sphere (SIS) the convergence 
profile as a function of radius R is: 

K(R) = ^ = -y(R), (ii) 



with i?Ein = 47r(<7sis/c) 2 -Dis/-Dos in radians and ctsis the 
lensing-inferred velocity dispersion which turns out to be 
very close to the stellar velocity dispersion of the lens 
galaxy (paper II). For a proper comparison, weak and 
strong lensing measurements have to be renormalized to 
the same source plane. Hence the Einstein radius given 
by strong lens modeling has to be rescaled by a factor 
^wl/wsl (see Table 1). Fig. 6 shows that after this 
scaling, but without fitting any free parameter, the strong 
lensing SIS models provide a reasonably good description 
of SLACS weak lensing data out to ~ 100 ft -1 kpc (with 
X 2 /dof ~ 22.8/9 ~ 2.5), and of the SDSS data beyond 
that (with x 2 /dof ~ 27.6/9 ~ 2.7). Two models are 
shown: one that neglects the non-linear relation between 
reduced shear and ellipticity (blue dashed) and one that 
takes this effect into account as well as the associated 
non-linear dependence on the source redshift distribution 
(solid green) . The two curves differ by less than the error 
bars of our measurements, showing that a linear relation 
between ellipticity and shear is a reasonable approxima- 
tion given the present statistical errors. This analysis 
shows that the total mass density profile of the SLACS 
lenses is very close to an isothermal sphere with velocity 
dispersion equal to the stellar velocity dispersion. Since 
the luminous component is steeper than isothermal out- 
side the effective radius, this finding implies the presence 
of an extended dark matter halo which is in turn shal- 
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Fig. 4. — Radial ellipticity profile in COSMOS left and SLACS right fields. Upper panels show the tangential component et 
profile and lower panels show the curl component e x . The signal is measured on both distortion-corrected stars (red triangles) 
and galaxies (black squares). For plotting convenience the innermost data point in the upper right panel (et in SLACS) falls 
outside the window. 



TABLE 2 

Measured excess surface density AS and shear 7 



Proj. Radius 


AS 


AS X 




It 


7x 




h kpc 


MVl /pc 2 


/iM Q /pc 2 


hM /pc 2 








3.3 


2307 


-1015 


1570 


0.339 


-0.054 


0.093 


5.8 


918 


-526 


614 


0.206 


-0.095 


0.078 


10.1 


115 


141 


264 


0.074 


0.010 


0.047 


17.6 


81 


-180 


153 


0.028 


-0.011 


0.029 


30.8 


232 


-36 


80 


0.064 


-0.005 


0.017 


53.9 


100 


36 


46 


0.019 


0.001 


0.010 


94.1 


90 


-45 


27 


0.016 


-0.010 


0.006 


164.5 


52 


52 


17 


0.010 


0.014 


0.005 


287.5 


60 


34 


17 


0.014 


0.008 


0.004 



The curl component AS X = 

lower than isothermal at similar radii 



it7x is statistically consistent with zero, cae is Iff error on both AS and AS> 

-ft and 7 X . 



ovy is la error on both 



5.1.2. Consistency with strong lensing and stellar 
dynamics 

A simple - although model dependent - way to com- 
pare on the same plot the mass measurement obtained 
with strong lensing, stellar dynamics and weak lensing 
is obtained in the following manner. For each radial bin 
we can define an effective weak lensing velocity disper- 



sion as the velocity dispersion of the singular isothermal 
sphere that reproduces the shear in that bin. The effec- 
tive weak lensing velocity dispersion profile is shown in 
Figure7 together with the average stellar velocity disper- 
sion determined from SDSS spectroscopy and the aver- 
age stellar velocity dispersion of the singular isothermal 
ellipsoid that best fits the strong lensing configuration. 

The figure illustrates the complementarity of the three 
mass tracers, stellar velocity dispersion well inside the 
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Fig. 6. — Radial shear profile around 22 SLACS strong lenses 
(black diamonds) expressed as the excess surface mass density 
AE(_R). SLACS bin are completely uncorrelated. The blue 
arrow deno tes t he mean effective radius of the lenses stellar 
component R c g ~ 5.3ft~ 1 kpc. The average critical density 
is E cr it — 4620 /iMq pc -2 . Previous results from S04 scaled 
to the same average square velocity dispersion are shown for 
comparison as red triangles (see § 4.4 for details). The curves 
show the shear profile expected from strong lensing SIS mod- 
els. This is not a fit. The solid green curve takes into ac- 
count non linear effects in the relation (e) = g = Ur)/{1 — wk) 
whereas the dashed blue line neglects non linearities. This 
shows that (e) = T07 is a good approximation for our pur- 
poses. 

and weak lensing in the context of a two-component mass 
model. 
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Fig. 5. — Mass reconstruction around stacked SLACS lenses 
in the same detector frame. The upper panel shows the con- 
vergence map (E-mode) whereas the lower panel shows the 
recovered convergence map after rotating galaxies by 45° (B- 
mode) which is consistent with a pure noise realization. 

Einstein Radius, strong lensing at the Einstein radius, 
and weak lensing outside the Einstein radius, as well 
as the dynamic range of the measurement, almost three 
decades in radius. The very close correspondence of the 
stellar and strong lensing measurement was discussed in 
paper II, and is confirmed here for a larger sample of 
lenses. This paper shows that, albeit with larger uncer- 
tainties, the weak lensing data show that the profile is ap- 
proximately flat for another two decades in radius. This 
is a qualitative statement as a full joint (strong+weak) 
lensing and dynamical analysis is needed to combine the 
three diagnostics properly. The three-pronged analysis is 
beyond the scope of this paper and is left for future work. 
In the rest of this paper we will focus on combining strong 
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Fig. 7. — Illustration of the complementarity of three mass 
tracers (stellar velocity dispersion, effective strong and weak 
lensing velocity dispersion) over almost three decades in ra- 
dius. Stellar velocity dispersion constrains the innermost re- 
gions (<C Reiu), strong lensing the region around the Einstein 
radius, while weak lensing is most effective beyond the Ein- 
stein radius. 
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5.2. Two component models 

In the rest of this section we perform a joint strong 
and weak lensing analysis of the data, in order to dis- 
entangle the mass profile of the stellar component and 
of the surrounding dark matter halo. For this purpose 
we adopt a simple two component model as detailed be- 
low. For simplicity, we assume that all lenses are at the 
center of their halo and none of them is an off-center 
satellite in a bigger halo. This approximation is well 
motivated by the galaxy-galaxy lensing results of M06 
who found that only a small fraction (< 15%) of massive 
ellipticals do not reside at the center of their host halo. 
Strictly speaking, the quantity measured by weak lensing 
is the projected galaxy-mass cross-correlation function 
rather than the actual shear profile of an individual halo. 
However, the interpretation of this cross-correlation func- 
tion within the successful framework of the "halo model" 
(e.g. Cooray 2006) allows one to disentangle the con- 
tribution of the proper halo attached to a given galaxy 
(1-halo central term), the halo of a more massive host 
galaxy (or group or cluster) if this galaxy is a satellite 
(1-halo satellite term), and the contribution due to clus- 
tering of neighboring halos about the main halo attached 
to the galaxy (2-halo term). However, the 2-halo terms 
only provide a significant contribution to the galaxy-mass 
cross-correlation function beyond a few Mpc (as com- 
pared to the outermost ~ 300 h~ x kpc radial bin probed 
here) and our lenses are massive ellipticals and thus most 
likely central galaxies. Therefore for the purpose of this 
analysis, and given the measured uncertainties, we can 
assume with little loss of accuracy that the measured 
shear profile is representative of the only surrounding 
main halo (see e.g. Fig. 1 of Mandelbaum et al. 2005). 

5.2.1. Model definition 

We model the stellar component with a de Vaucouleurs 
density profile (de Vaucouleurs 1948; Maoz & Rix 1993; 
Keeton 2001) The effective radius of the stellar com- 
ponent is fixed to the ACS surface photometry. Thus 
the only free parameter needed to measure the luminous 
component is the average stellar mass-to-light ratio Ty. 
The dark matter halo is assumed to be of the NFW form 
(Navarro et al. 1997; Bartelmann 1996; Wright & Brain- 
erd 2000) in which the density reads 



p(r) = p c 6 c 



1 



(12) 



with r s the scale radius and p c the Universe critical den- 
sity. The term 5 C = f c 3 / [ln(l + c) - c/(l + c)] relates 
the so-called concentration parameter c = r&/r s . We as- 
sume an overdensity of A = 119 so that r& can be con- 
sidered the "virial" radius (Bryan & Norman 1998). This 
definition agrees with those of Hoekstra et al. (2005) and 
Heymans et al. (2006a) for our assumed ACDM cosmol- 
ogy. Because statistical errors are still large with only 22 
SLACS lenses used here, we lack the sensitivity to con- 
strain the dark matter profile in detail. Thus, instead 
of fitting a free concentration parameter, we assume the 
general relation observed in numerical simulations: 



M v 



1 + z V 8.12 x lO 12 /^ 1 ]^ 



-0.14 



(13) 



(Bullock et al. 2001; Eke et al. 2001; Hoekstra et al. 
2005). In addition, we are not able to constrain the virial 
mass of each lens individually, so we need to assume a 
scaling relation between virial mass and V band luminos- 
ity of the form: M v ; r = TyLy. Note that we check that 
assuming a steeper relation M v ; r oc L v 5 (Guzik & Sel- 
jak 2002; Hoekstra et al. 2005; Mandelbaum et al. 2006) 
yields comparable results because the SLACS lenses span 
a narrow range in luminosities, with 0.2 dex rms around 
(Ly) — 5.70 x 10 10 /i~ 2 L Q . In conclusion, our model 
has only two free parameters, the mass-to-light ratio of 
the luminous component Ty and the virial mass-to- light 
rato Ty. 

Having defined the model, we now define the merit 
function that will be used to determine the best fit- 
ting parameters with their uncertainties. Detailed strong 
lensing analysis of multiply imaged sources have shown 
that extremely tight constraints can be set on the Ein- 
stein radius of individual lenses (e.g. paper III), with typ- 
ically a few percents relative accuracy o\R Ein /i?Ein — 5%. 
This can be interpreted as a surface mass measurement 
since the mean density S within this radius is by defi- 
nition equal to the critical density S C rit- Therefore for 
each lens we are able to write: 



£(< i?Ein) = £ crit (2i, z s ) ■ 



(14) 



The relative error on i?Ein translates into cr^/T. ~ 5%. 
We can thus define a strong lensing merit function: 



Xsl 




^>*,i(Yy) — T,DM,i(Ty) 



(15) 



where subscripts * and DM stand for luminous and dark 
matter components, evaluated at position i?Eini- 

In a complete analogy, we define the weak lensing merit 
function: 

N rbin 



Xwl ^ ^ 



1 



J=1 AEj 



AS, 



lens 



J2 (AE.,y (T V ) + AS DM ,y(7y)) 



i=l 



(16) 

where iV r bi n is the number of radial bins rj at which 
ASj = AS(r-j) is obtained from the weak lensing data 
shown in Fig. 6. 

In the next section we will derive the best fitting 
{Ty, Ty} values that minimize the total \ 2 = Xsi + Xwl- 

5.2.2. Results 

The upper panel of Fig. 8 shows the radial pro- 
file of the shear for the best fit model together with 
weak lensing data points. The fit is excellent with a 
X 2 /dof = 29.1/31 ~ 0.94. We see the detail of the con- 
tribution of stellar and dark matter components. This 
joint strong+weak lensing analysis allows to disentangle 
the contribution of each. Because the latter component 
is less concentrated and more extended than the former, 
there is a radial range R ~ 20 h^ 1 kpc at which surface 
mass density flattens. This implies a fast drop in the 
shear profile AS (J?) at that scale which is quite easy to 
detect. 
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Fig. 8. — Upper panel: Shear profile (i.e. AE) for the best 
DM + de Vaucouleurs profile. The contribution of each mass 
component is detailed (green and blue for stars and DM re- 
spectively). The thickness of the red curve codes for the la 
uncertainty around the total shear profile. Uncertainties are 
very small below 10 kpc because of strong lensing data that 
cannot be shown here. Lower panel: Similar coding for the 
three-dimensional density profile p(R). The transition be- 
tween star- and DM-dominated mass profile occurs close to 
the mean effective radius (blue arrow). The total density 
profile is close to isothermal over ~ 2 decades in radius. 

In the lower panel of Fig. 8 we show the correspond- 
ing three-dimensional mass density profile p(R). This 
profile is close to isothermal although it is made of two 
components which are not isothermal. The components 
combine to make an almost isothermal density profile 
at scales R < 100 ft. -1 kpc with a transition from star- 
dominated to dark matter dominated profiles occurring 
close to the effective radius. 

The best fit NFW + de Vaucouleurs lens model yields a 
stellar mass-to- light ratio Ty = 4.48 ±0.46 /iM Q /L Q and 
virial mass-to-light ratio Ty = 246^g7 1 HMq/Lq. This 
translates into a virial to stellar mass ratio M v ; r /M* = 
54^21- Note that stellar mass-to-light ratio depends in 
the same way on h as virial mass as they are inferred 
from lens modeling and not from stellar evolution models. 
Thus M v i r /M* is independent of h. Fig. 9 shows the 
68.3, 95.4 and 99.7% confidence levels contours for the 
best fit model parameters. A model having a constant 
mass-to-light ratio (i.e. M vil /Ly = 0) is ruled out at 

- 4(7. 




Fig. 9. — Confidence levels around model parameters (1, 2 
and 3 a contours) for the relation between total virial mass- 
to- light ratio Ty ~ M v i r /Lv and stellar mass-to-light Ty = 
M t /Ly ratio. Given the sample mean luminosity (Ly) = 
5.7 x 10 10 Ii~ 2 Lq, we find a mean sample stellar mass (M*) = 
2.55 ± 0.27 x 10 11 h _1 M© and virial mass (M vlr ) 
10 la A _1 M B . 
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Given the sample mean luminosity (Ly) = 5.7 
10 10 h~ 2 L Q , we find a mean sample stellar mass (M*) 
2.55 ± 0.26 x 10 11 Ii-^r and virial 
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(M vir ) = 

10 12 ft, _1 M . This translates into a mean virial 
(resp. scale) radius r vir = 3931*3 kpc (resp. r 8 — 
58 ± 8/i _1 kpc). We note that the virial radius is typ- 
ically larger that our field of view, and therefore virial 
masses rely on extrapolations of our results. Therefore, 
we also present more robust measurements like the pro- 
jected and three-dimensional mass within a reference ra- 
dius R — 200 /i -1 kpc. Lens modeling yields M 3 d(< 
200 h' 1 kpc) = (8.1 ± 1.8) x 10 12 h^M© and a projected 
mass M 2D (< 200 h' 1 kpc) = (10.8 ± 2.7) x 10 12 /i^Mq 
(68% CL errors). 

To compare with local measurements we convert our 
Ty mass-to-light ratio to the rest-frame B band. As- 
suming a typical (B — V) — 0.96 color for Ellipticals 
(Fukugita et al. 1995), (B - V) = 0.65 and a Hubble 
constant h = 0.7, one finds T B = 4.17 ± O.44M /L . 
Using paper II, Treu & Koopmans (2004) and similar 
findings (Treu & Koopmans 2002; Treu et al. 2005; van 
der Wei et al. 2005; di Serego Alighieri et al. 2005) for 
the passive evolution of massive early-type galaxies, we 
get a redshift zero B band stellar mass-to-light ratio 
~ 5.81 ± 0.61 which is statistically consistent with lo- 
cal estimates such as T b = 7.8 ± 2.7 and = 7.1 ± 2.8 
from Gerhard et al. (2001) and Trujillo et al. (2004) re- 
spectively. The low value of M*/L is also in broad agree- 
ment with stellar evolution models of Bruzual & Chariot 
(2003), although detailed comparisons depend on the as- 
sumed Initial Mass Function (IMF). 

Our modeling can be used relate the V band luminosity 
within the Einstein radius Ly(< i?Ein) and the fraction 
of dark matter in the same projected radius /i)M,2d(< 
-REin)- Using Eq. (14), for a given stellar mass-to-light 
ratio each lens must verify 

/DM,2d(< -REin) = 1 - YyLy(< #Ein) /Scrit • (17) 

Fig. 10 shows the inferred projected /dm. 2d using our 
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best fit Ty = 4.48 /i _1 M Q /L0. The mean dark matter 
fraction within the Einstein radius (/dm, 2d) = 0.37 ±0.04 
with 18% rms scatter. Extrapolating to the effective ra- 
dius, about half of the projected mass is in the form 
of dark matter. The result from the NFW + de Vau- 
couleurs parameterization is shown as the solid curve 
which matches the data points well (see also papers II 
and III). This parameterization also allows the depro- 
jected dark matter fraction to be calculated, and is found 
to be ~ 27% (dotted blue line) within R c g. The local de- 
projected DM and stellar densities are of the same order 
at that radius. 




1 10 
radius [kpc/h] 

Fig. 10. — Projected dark matter fraction in SLACS lenses 
as measured from strong lensing at the Einstein radius (data 
points) and our constraints on the stellar M*/L ratio. The 
Einstein radius is expressed in physical units and the solid 
curve shows the best fit profile of the dark matter fraction 
/dm, 2d inferred from the parametric NFW + de Vaucouleurs 
modeling of strong+weak lensing data. The mean effective 
radius is shown as a blue arrow. We see that /dm, 2d (< -Reft ) — 
50% of the projected enclosed mass is in the form of dark 
matter within the effective radius. 

It may seem that the two data points with /dm. 2d — 
are responsible for our inferred low M»/L. Since these 
two lenses have the most elongated stellar component 12 , 
the assumption of circular symmetry may break down for 
them. Furthermore, if they have a disk component with 
younger starsq, this could reduce their global M*/L with 
respect to that of pure spheroidal systems. However, re- 
doing the strong+weak lens modeling without these does 
not change the stellar mass-to-light ratio to significantly 
(M*/Ly =4.90 ±0.53 ft- ^q/Lq). 

5.3. Comparison to previous galaxy- galaxy lensing works 

We first compare our findings to the SDSS weak 
galaxy-galaxy lensing analysis of M06 who defined M cen t 
as the mass enclosed in a sphere within which the mean 
density is 180 times the background density, similar to 
our definition. Our lens sample should be compared 
to the sm7 stellar mass bin for early-type galaxies with 
(M*) = 39.6 x 10 10 M Q . It would also lie between 16b 
and 16f early-type galaxy luminosity bins with a typ- 
ical (L) ~ 6.3L*. To convert the rest frame V band 
to the SDSS r' band we use V - r' ~ 0.36 for Ellip- 
ticals (Fukugita et al. 1995). Since at the exponen- 
tial tail of the luminosity and mass functions one is 

12 g* = (fe/a)* ~ 0.5 whereas the other lenses have a mean 
g» = 0.81 and dispersion 0.08 



highly sensitive to the scatter in the mass-luminosity 
relation, M06 applied corrections calibrated into simu- 
lations (Mandelbaum et al. 2005) whereas our analy- 
sis does not attempt to correct for this effect. For the 
16f bin, M ccnt /L r = 674± 2 ,J 3 I /iM Q /L as found by M06 
yields M ccnt /L v = 896 ± 275 /iM Q /L Q (95% CL). Like- 
wise the sm7 bin of M06 yields M cent /M, = 256±gg (95% 
CL and using h — 0.7). These values are significantly 
larger than our results. However if we apply a simi- 
lar correction as these authors our virial mass should 
be raised by 66%, yielding M vil /L v = 408±^ hM Q /L Q 

and M v i r /M* = 89+35- This latter correction thus brings 
our findings into statistical agreement. The small differ- 
ence might be due to our inability to probe the very outer 
parts of halos and efficiently constrain virial masses. 
However, we emphasize that the availability of strong 
lensing constraints put tight constraints on the column 
density enclosed by i?Ein, which means that the outer 
parts of halos cannot contribute much in making lens 
galaxies critical, or perhaps to the fact that we fit for 
M*/L, while M06 use the value determined from stel- 
lar population synthesis models. Issues related to com- 
parisons between Guzik & Seljak (2002) and M06 re- 
sults are addressed in the latter paper. In any case, we 
find a much better agreement between our M v ; r mea- 
surement and those of Guzik & Seljak (2002) which give 
M v i T /Lv — 296 ± 51ft-M Q /L0 after matching our lens 
sample selection. 

We now compare our analysis with the Hoekstra et al. 
(2005) results. Since we chose to match their definition 
of virial mass and concentration, we expect comparisons 
to be easier although the mean redshift of their lens sam- 
ple is ~ 0.32. Our lens sample would passively brighten 
to (L v ) ~ 6.3 x 10 10 h- 2 L at z = 0.32 which is - 2.45 
times brighter than their higher luminosity bin having 
Ly ~ 2.45 x 10 10 /i~ 2 Lq. Therefore we need to ex- 
trapolate their results using their M w -„ cx L 15 relation. 
They find M vir /Ly ~ 253±|| /iM Q /L . This value is 
statistically consistent with ours. However, a detailed 
comparison is made difficult due to the fact that the au- 
thors mix early- and late-type galaxies and they avoid 
lens galaxies in dense environments. Using stellar evo- 
lution models, they estimate the virial to stellar mass 
in their reddest subsample (B — V) rcs t ~ 0.95 to be 
M vir /M* ~ 43 ± 6 for a scaled Salpeter IMF (Bell & de 
Jong 2001) or M vir /M* ~ 27 ±4 for a PEGASE Salpeter 
IMF. The Scaled Salpeter IMF hypothesis turns out to 
be in better agreement with our measurements. In ad- 
dition it predicts stellar mass-to-light ratios Tg ~ 4.5 
(Solar units) closer to our estimates than the PEGASE 
IMF which predicts ~ 6.5. 

Finally, Heymans et al. (2006a) measure virial masses 
of lens galaxies in the range 0.2 < z\ < 0.8 in a narrow 
range of luminosity L r — 2 Ax 10 10 h~ 2 L© with a 0.2 dex 
dispersion about this mean. This sample is dominated by 
early-type galaxies. Again, extrapolation to our sample 
mean luminosity is somewhat uncertain but, using the 
M v i r cx L v 5 scaling, their results give M v i r /M* = 76 ± 
25 which is in excellent agreement with our M v i r /M* = 
54tg. 

These comparisons show that our virial mass estimates 
are in good agreement with other studies after extrap- 
olation of our constraints on the radial shear profile 
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(ending around ~ 300/i _1 kpc) out to the virial radius 
~ 480ft.~ 1 kpc. Comparing with other results obtained 
for less massive systems on average increases uncertain- 
ties. However our results on the halo virial masses are 
well consistent with this ensemble of results above as they 
lie in between them. In addition, the measured shear 
profile remarkably match those of S04 and M06 in the 
radial range 30 < R < 300/i _1 kpc. This gives a valu- 
able support to the validity of our results and the control 
of systematic effects. 

6. DISCUSSION 

Our joint weak+strong lensing modeling of SLACS 
lenses with a two-component mass model allowed us to 
successfully disentangle the contribution of each giving 
sensible results for both the stellar mass-to-light ratio 
M*/L v = 4.48 ± 0.46 /iM©/L© and virial mass-to-light 
ratio M v i r /Ly = 246tg° 1 /iM©/L©, in good agreement 
with other studies. Assuming a NFW and de Vau- 
couleurs forms for each density profiles provides a good 
description of the data (x 2 — 0.94 per degree of freedom). 

This analysis shows that the total density profile is 
close to isothermal out to ~ 100 h^ 1 kpc. It is now a well 
established from SLACS (paper III) and earlier strong 
lensing studies (e.g. Treu & Koopmans 2002; Rusin et al. 
2003; Rusin & Kochanek 2005; Treu & Koopmans 2004) 
that the total mass profile of lens galaxies is close to 
isothermal (p oc r~ 2 ) within R e g. In paper III this result 
was established by combining strong lensing and stellar 
kinematics. The present analysis extends and strength- 
ens this results as we find that the dark matter and stellar 
components combine themselves to form an isothermal 
total density profile well beyond the effective radius. 

We find that the transition from a star-dominated to a 
dark-matter-dominated density profile must occur close 
to R e s. This peculiar transition is also observed by Treu 
& Koopmans (2004) which combined strong lensing and 
stellar kinematics in higher redshift lenses. Similar re- 
sults can be found in Mamon & Lokas (2005a, b). In ad- 
dition, the "mean" fraction of DM fdm,SD(< -Roff) ~ 30% 
is in excellent agreement with local estimates (Kronawit- 
ter et al. 2000; Gerhard et al. 2001; Borriello et al. 2003), 
and more recently from the SAURON project (Cappel- 
lari et al. 2006). At this point, it is noteworthy to note an 
important result of paper II, that is, strong lensing galax- 
ies have similar internal properties as normal early-type 
galaxies in terms of their location in the Fundamental 
Plane. Our findings can thus be generalized. 

We emphasize that we did not investigate other pa- 
rameterizations for the DM halo. For instance, a steeper 
profile (pdm oc r~ a with a > 1) could possibly arise 
from the adiabatic contraction (Blumenthal et al. 1986; 
Gncdin et al. 2004) of a NFW halo which is found in pure 
dark matter N-body cosmological simulations made with- 
out taking into account the complex physics of baryons. 
With the small sample of lenses we are considering here, 
we are not sensitive yet to that level of detail in the 
inner slope of the assumed DM profiles. We note how- 
ever that the rather low values of M^jLy we find would 
make it unlikely to have a dark matter halo much steeper 
than NFW (p oc r _1 at the center) (see also Borriello 
et al. 2003; Humphrey et al. 2006). Since it is reason- 
able to assume that baryons somehow perturb the DM 
halo within the effective radius, we emphasize that our 



successful NFW parameterization should rather be con- 
sidered as a fitting formula for the perturbed halo. In 
future papers, with the complete observed lens sample 
and spatially resolved measurements of stellar kinemat- 
ics, we plan to determine with unprecedented accuracy 
the inner slope of the DM profile below R a g and at the 
same time M^jLy . 

The inner regions of lens galaxies can be considered 
as representative of the whole parent sample of early- 
type galaxies as shown in paper II. Although there is no 
firm observational evidence, it is thought that environ- 
mental effects may bias the population of lens galaxies 
since extra convergence coming from surrounding large 
scale structure may boost lens efficiency while leaving 
the internal dynamics of lens galaxies unchanged (see 
also Keeton & Zabludoff 2004; Fassnacht et al. 2006). 
In the present analysis we address the issue of whether 
lenses are representative of the overall population at 
larger radii, by comparing our weak lensing results with 
those obtained for non-lens samples. The present anal- 
ysis shows that, on intermediate scales R e s < R < 
300/i _1 kpc, SLACS massive lenses have the same shear 
properties as normal Ellipticals as found by M06 or S04. 
We find a good agreement between our virial mass es- 
timates and semi-analytic predictions like those devel- 
oped in the "halo model" (Mandelbaum et al. 2005, 2006; 
Cooray 2006). Any systematic environmental effect able 
to perturb shear measurement on those scale is below 
our present observational uncertainties. When the ACS 
follow-up is finished, the weak lensing analysis will pro- 
vide important new information on the environment of 
strong lens galaxies. This issue is deferred for a future 
work. 

In terms of the internal structure of early- type galaxies, 
the present analysis strengthens and extends the results 
presented in paper III and gives additional support to the 
picture proposed there for their formation. The isother- 
mal density profile must be produced at early stages of 
their evolution process (z > 2) by merging/ accretion in- 
volving dissipative gas physics to quickly increase the 
central phase space density since pure collisionless sys- 
tems would rather develop shallower density slopes p oc 
rl" 1 - 1 - 5 ! (e.g. Navarro et al. 1997; Moore et al. 1998; 
Ghigna et al. 1998; Jing 2000; Navarro et al. 2004). Once 
isothermality is set, early- type galaxies may evolve pas- 
sively, or may keep growing quiescently via collisionless 
"dry" mergers or by accretion of satellites, which pre- 
serve the inner density profile of collisionless materials 
(stars and/or DM). See discussion in paper III and ref- 
erences therein for further details. 

7. SUMMARY & CONCLUSION 

We demonstrate that, using deep ACS images, it is pos- 
sible to measure a weak lensing signal for a modest sam- 
ple of massive early- type galaxies (a v ~ 250 km s -1 or 
logM*/M© ~ 11.35): with only 22 lenses we are able to 
detect shear signal at 5cr significance. Key to this success 
is the large density of well-resolved background sources 
afforded by ACS, which beats down shot noise and re- 
duces the problem of dilution by lens satellites since back- 
ground sources greatly outnumber satellites. In addition, 
these lenses are very distant from one another and hence 
completely statistically independent. Furthermore, spe- 
cial care has been taken to control systematic errors, us- 
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ing the most advanced techniques to model and correct 
for the ACS PSF and other instrumental effects. By an- 
alyzing a sample of 100 blank fields from the COSMOS 
survey we show that residual systematics in the shear 
measurement are less than 0.3%. 

Although weak lensing alone can provide interesting 
results on lens density profiles at intermediate scale, the 
great power and originality of this work is the combi- 
nation of strong and weak lensing constraints. Mod- 
eling weak and strong lensing in massive ((c^) 1//2 ~ 
248kms _1 ) SLACS galaxies as a sum of a stellar (de 
Vaucouleurs) plus dark matter halo (NFW) components, 
we could disentangle to contribution of the two compo- 
nents overall mass budget. The main results of this joint 
analysis can be summarized as follows: 

• The total density profile is close to isothermal from 
0.5i?cff ;$ R < 100-Rcff although neither the stellar 
nor the dark matter density profile is isothermal. 

• The transition from star- to dark matter- 
dominated density occurs at ~ i? ff, leading to a 
dark matter fraction within this radius fdm3D(< 
i? off ) = 27 ± 4%. 

• The best fit stellar mass-to-light ratio is M^/Ly = 
4.48 ±0.46 HMq/Lq in agreement with local results 
and old stellar populations. 

• The best fit virial mass-to-light ratio is M v ; r /iy = 
246^gy 1 ftM©/L Q in agreement with galaxy-galaxy 
weak lensing results of non-lens galaxies. We found 
a mean virial mass and radius (M v ; r ) = 14^5 x 
HF/T^M© and r vir = 393 ± SO/i^kpc respec- 
tively. 

• The agreement with other weak lensing studies 
shows that the outer halos of lenses and non-lenses 
are consistent within the errors. In other words, 
if lens early-type galaxies live in peculiar environ- 
ments, their effect on the shear profile down to 
300 h^ 1 kpc from the lens center are below our sta- 
tistical errors. 

We forecast a ~ 10er detection by the end of the ongo- 
ing deep follow-up imaging with HST/ACS. When com- 
pleted, the SLACS sample of lenses with well resolved 
kinematics will provide valuable constraints on stellar 
populations and density profiles of both stellar and dark 
matter components down to several hundred kiloparsecs 
of the lens center, thus allowing to address internal prop- 
erties of lens galaxies as well as the effect of their envi- 
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ronmcnt. To complete the picture on early- type galax- 
iess and structure formation, it is important to extend 
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